Optical investigations of quantum-dot spin dynamics 
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We have performed all-optical measurements of spin relaxation in single self-assembled InAs/GaAs 
quantum dots (QD) as a function of static external electric and magnetic fields. To study QD spin 
dynamics we measure the degree of resonant absorption which results from a competition between 
optical spin pumping induced by the resonant laser field and spin relaxation induced by reservoirs. 
Fundamental interactions that determine spin dynamics in QDs are hyperfine coupling to QD nu- 
clear spin ensembles, spin-phonon coupling and exchange- type interactions with a nearby Fermi sea 
of electrons. We show that the strength of spin relaxation generated by the three fundamental 
interactions can be changed by up to five orders of magnitude upon varying the applied electric and 
magnetic fields. We find that the strength of optical spin pumping that we use to study the spin 
relaxation is determined predominantly by hyperfine-induced mixing of single-electron spin states 
at low magnetic fields and heavy-light hole mixing at high magnetic fields. Our measurements allow 
us to determine the rms value of the hyperfine (Overhauser) field to be ~15 mTesla with an electron 
g-factor of ge=0.6 and a hole mixing strength of \ef = 5 x 10~*. 

PACS numbers: 



I. INTRODUCTION 



A single quantum dot (QD) electron spin is a fun- 
damental physical system which allows for a controlled 
study of confined spin dynamics in the solid-state. In 
contrast to higher-dimensional semiconductor structures 
QD spins have been demonstrated to posses long relax- 
ation and coherence times exceeding 20 msec and 10 fj,sec 
respectively. These findings along with demonstration of 
single spin manipulation and read-out have strengthened 
the proposals for using QD spins as physical representa- 
tion of qubits in quantum information processing 
The prolongation of spin relaxation times for QD spins 
stems from a drastic reduction in spin-phonon coupling 
mediated by a combination of electron-phonon and spin- 
orbit interactions and suppressed by strong quantum con- 
finement of electrons. As a consequence additional spin- 
reservoir interactions such as hyperfine coupling to QD 
nuclear spins and exchange-type (co-tunneling) coupling 
to a nearby Fermi-sea become prominent in determining 
the spin dynamics in QDs. 

Here we study the dynamics of an electron spin con- 
fined in a self- assembled InAs/GaAs QD which is in 
turn embedded in a Schottky heterostructure. In order 
to assess the relative importance and external field de- 
pendence of the three elementary spin-relaxation mech- 
anisms we use the degree of resonant absorption as a 
measure: since the degree of absorption i.e. strength of 
the trion-resonant light scattering is determined by com- 
peting optical spin pumping (OSP) [1] and spin relax- 
ation induced by spin-reservoir interactions, the strength 
of spin relaxation can be inferred from the absorption 
measurement. First we demonstrate that at low mag- 
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netic fields (up to 1 Tesla) spontaneous spin-flip Ra- 
man scattering that allows for one-way pumping into 
the optically-uncoupled spin-state is predominantly me- 
diated by a mixing between the electronic spin states in- 
duced by the fluctuating hyperfine nuclear (Overhauser) 
field. Next, we show that upon varying the external gate 
voltage by about 50 mV the spin relaxation due to ex- 
change coupling to the nearby Fermi-sea of electrons can 
be changed by as much as five orders of magnitude. Fi- 
nally, we show that at the high magnetic field regime 
(1 to 10 Tesla) spin pumping is due to heavy- light hole 
mixing and spin relaxation is dominated by phonons in 
conjunction with spin-orbit interaction. 

Before proceeding we note that major advances in un- 
derstanding relaxation and decoherence of single confined 
electron spins have already been achieved in electrically 
defined QDs. By implementing a single-shot electrical 
read-out of a QD spin Elzerman et al have shown that 
spin lifetimes in electrically defined QDs can reach up to 
^1 msec even at elevated magnetic fields of 8 Tesla @ 
and very recent measurements have revealed a relaxation 
time of 170ms at 1.75 Tesla |t6J]. Similarly in double QDs 
Johnson et al have investigated hyperfine-induced triplet- 
singlet relaxation ^ Petta et al have demonstrated co- 
herent manipulation of singlet-triplet states [8] and Kop- 
pens et al have shown detection and control of hyperfine- 
induced singlet-triplet mixing Q and Rabi oscillations 
using microwave pulses (lo| . Further measurements on 
InAs/GaAs self-assembled QD ensembles have revealed 
Ti times exceeding 20 msec at a magnetic field of 4 Tesla 
and a temperature of 1 Kelvin [llj . 

This paper is organized as follows: In section [TT] we in- 
troduce the coupling of the localized spin to nuclear spins, 
charge reservoir and phonons. Section [1111 then theoreti- 
cally describes the QD spin dynamics in the framework of 
the trion four-level system with spin-reservoir coupling. 
In section IIVI we present our experimental results ob- 
tained with single QD absorption spectroscopy in distinct 



2 



regimes of external electric and magnetic fields where dif- 
ferent interactions dominate. Finally section IVl gives an 
overview on the above-mentioned interactions together in 
a self-contained picture before the conclusions in section 
IVII For information about our sample structure and ex- 
perimental techniques we direct the reader to Appendix 

m 

II. INTERACTIONS OF A SINGLE CONFINED 
SPIN 

A. Nuclear Spins 

The interaction of a localized electron spin with a sur- 
rounding nuclear spin ensemble can be written in the 
form of the Fermi contact interaction which yields 

i 

The sum runs over all nuclei i in the lattice, vq is 
the volume of an InAs unit cell ij^iRi) the electron 
envelope wavefunction at the ith nucleus and li and 
a are the spin operators of nuclear and electron spin. 
Ai = {2iJLQgoiJLBHi/'liIi) \uc{Ri)\^ is the hyperfine coupling 
strength and reflects the electron density described by 
the electron Bloch wavefunction Uc{Ri) at the site of the 
nuclei, /is is the Bohr magneton and /i^ the nuclear mag- 
netic moment /io is the permeability of vacuum and 
the free-electron g-factor. 

In order to estimate the total number of nuclei within 
the spread of the electron wavefunction we use the di- 
mensions of the QD: InAs and GaAs have Zincblende- 
type lattice with a lattice constant of 6.06 A and 5.65 A. 
There are four Arsenic (nuclear spin Ias = 3/2) and four 
Indium (Iin=9/2) or Gallium (lGa=3/2) atoms in a fee 
unit cube. Taking this into account and assuming that 
the QD creates a box-like confinement that equals the 
dimensions of the QD 20nmx20nmx5nm) the num- 
ber of As atoms is on the order of 4x10'* and one has to 
add the same number of In or Ga atoms depending on 
the composition of the QD such that the total number of 
nuclei interacting with the QD spin can be taken to be 
N = lO'' to 10^ 

We note that the Fermi contact interaction relies on a 
finite value of the Bloch wavefunction at the sites of the 
nuclei. Due to the s-like symmetry of their wavefunction 
electrons are susceptible to this interaction whereas the 
p-symmetric holes are not. For holes it is only possible 
to interact with nuclear spins via the much less efficient 
dipole-dipole interaction. 

The hyperfine interaction constitutes a special case 
among the three spin-reservoir interactions discussed in 
this work as our experiments suggest that its dominant 
effect is not described by a textbook system-reservoir in- 
teraction in the Born-Markov approximation as it is the 
case for phonon and exchange coupling. In the optical 



measurements presented here nuclear spins mainly act by 
exerting a quasi-static magnetic field ( Overhauser field) 
with rms- value on the order of i3nuc=15 mTesla. This 
field leads to a Rabi-type slowly varying coherent mixing 
of the spin ground states which can also be understood 
as precession of the electron spin in the nuclear magnetic 
field; in this context slowly varying means that the cor- 
relation time of the hyperfine field fluctuations (~lms) 
is much longer than the precession time (~lns) of the 
electron spin in this field. 

In contrast to the ground states the excited states re- 
main unchanged as the hole is not susceptible to the nu- 
clear magnetic field and the two electrons form a singlet 
which is immune to magnetic field variations too. The 
mixing of the ground states if strong enough would thus 
lead to a fluctuating observable splitting of the excitonic 
transitions however it turns out that its magnitude is less 
than the broadening of these transitions due to their ra- 
diative lifetime and therefore cannot be resolved in laser 
scanning absorption measurements even with experimen- 
tal resolution much better than the transition linewidth. 

In order to understand the effect of the hyperflne field 
in optical experiments we consider two regimes: First an 
external magnetic field with strength smaller than the hy- 
perfine field is applied or the external field is completely 
absent. Hence the direction of the total magnetic field 
seen by the electron spin is fully random after a nuclear 
field correlation time. As we will see in section IIIII fast 
bidirectional OSP will be the consequence inducing effi- 
cient spin relaxation dominating over other mechanisms 
as the one predicted in Ref. [l^ . In the second regime the 
applied external field is much stronger than the hyper- 
fine field. In this case the electron spin mainly sees the 
external magnetic field along the z-axis and the hyper- 
fine field only leads to small fluctuations of the nuclear 
field vector. In this regime the light-induced spin relax- 
ation is slow and other mechanisms such as phonons are 
dominant. 

We will now analyze the effective hyperfine field in de- 
tail. In the absence of dynamical nuclear spin polariza- 
tion (DNSP) scenarios the action of nuclear spins upon 
the localized spin dominates the reverse action due to a 
much larger number of degrees of freedom on the side of 
the nuclear spins ensemble. Therefore the Hamiltonian 
(P) is reduced to an effective magnetic field seen by the 
QD spin which is commonly referred to as Overhauser 
field 



B^-i^C^i.) (2) 

where A is an average spin-nuclei coupling constant and 
ge is the QD electron g-factor. 

For simplicity we will in the following treat the hyper- 
fine field as a purely classical field Bjv {t) with correlation 
time Tcorr '^-^Ims. The correlation time is expected to be 
similar to the decay time of nuclear spin polarization in 
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the presence of a QD electron and absence of external 
magnetic field, as measured in i?nuc refers to the 

rms- value of the Gaussian distribution as defined by 



and we define 



9ef^BBN.xy{t) 



(12) 



f{BM) 

which yields 



1 



i?Lc(27r)3/2 



exp(-|B;v|V2i?Lc) 



(Bjv(t)) = 
(|B^(t)P) = 3B, 



(3) 



(4) 



2 

nuc 



Here ( ) denotes the time average over many correlation 
times. 

As our Bjv(i) is quasi-classical, we treat the BN,i{t) 
with i = x,y,z as independent random variables each 
one following a Gaussian distribution function given by 



/27r 



exp 



(5) 



and /(Bat) is then given by the product of the Gaussian 
distributions of its three spatial components 



/(Sjv) = n /(^^■') 



Clearly 
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can be written in the form (similar to |12l. Il5|) 
bo 



B,, 



(6) 



(7) 



(8) 



with bo a parameter characterized by the species of nuclei 
and the composition of the QD and N the number 
of nuclear spins interacting with the QD spin. 

The QD com position is taken to be 90% InAs and 10% 
GaAs yielding I{I +1) = 13.2 when averaging over the 
different nuclear species [l6|. Similarly we obtain = 
2500/ieV2 which yields bo = 3.0 Tesla. Using © with 
N^IO" to 10^ nuclear spins we obtain for our QDs 



B„ 



9.5 mT to 30 niT 



We can now rewrite (H]) as 



(9) 



(10) 



Here the component of the nuclear field along the z-axis 
BN.z{t) only leads to Zeeman splitting whereas the in- 
plane components induce a mixing of the | "f) and | J,) 
states. The in-plane hyperfine field is 



We note here that our measurement time (typically 10 to 
100 ms) is longer than the correlation time of the nuclear 
field; i.e. for each measured data point we expect that we 
average over many configurations of the nuclear magnetic 
field El- 



B. Coupling to electron spin reservoir 

Another interaction mechanism arises from the pres- 
ence of the Fermi sea in the back contact (rcf. to Ap- 
pendix [X]) that couples to the QD via the tunneling 
barrier. It is well known that exchange interaction of 
a confined spin with an electron spin reservoir leads to 
co-tunneling [l7| [s^ and at temperatures lower than the 
Kondo temperature Tk to the formation of a Kondo sin- 
glet p,iia[20,[2i|. 

This interaction can be written as 



charge — / ^ 
k,k' 



%t,fe(eLi^QD,TeQD,iefe',T + cc.) (13) 



where 



QD,a 



and eQD,cr are the creation and destruc- 
tion operators for an electron with spin a in the QD and 
similarly in the reservoir, gt^k is the tunneling matrix 
element, which is linked to the tunneling rate Ftunnoi by 
Fermi's Golden rule Ftunnei = ^\gt,k\'^ p{E) with p{E) 
being the density of states in the back contact. 

Fig. [T^a) shows the energies of the empty, singly and 
doubly charged QD state as a function of gate voltage 
[22] . Which state has lowest energy obviously depends 
on the gate voltage and the QD attempts to reach it by 
either attracting or repelling electrons from or into the 
reservoir. Clearly there is a range of voltages {single elec- 
tron charging plateau) where it is energetically favorable 
for the QD to accomodate a single electron, marked by 
the shaded region in the figure. At the points A and B 
two charging levels are degenerate and fast exchange of 
the QD electron with the reservoir can take place, only 
limited by the tunneling rate. The real gate voltages Va 
and Vb that need to be applied in order to reach points 
A and B can vary from dot to dot depending on its con- 
finement properties. 

In other words the QD is singly-charged for Va < Vg < 
Vb with Vg the gate voltage. We define the plateau center 
14 



Vb-Va 



(14) 



The gate voltage detuning is 



Bli^xyit) = Bjf.^{t) + B%y{t) 



(11) 



AVg =Vg~V 



(15) 
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FIG. 1: (Color online) (a) Stability diagram of the QD ground states neglecting spin; Energies of the zero one and two-electron 
QD as a function of gate voltage. Crossover points are marked A and B. E12 denotes the gate volt age- dependent energy 
difference between the singly-charged and the doubly-charged state or charging energy, (b to d) The QD can exchange its 
single electron with the charge reservoir via a virtual empty or two-electron (shown here) state. When one of the two singlet- 
electrons tunnels out it leaves the remaining QD spin in a mixed state equivalent to spin relaxation. Ftunnci marks the tunneling 
rate through the 35-nm GaAs barrier e the detuning from the Fermi energy ef- Ez is the Zeeman splitting and E12 the energy 
required to charge a second electron. 



The schematic co-tunneling process is depicted in Fig. 
[IJb-d). The initial state is characterized by a QD with 
a single spin-down electron and Coulomb blockade pro- 
hibits tunneling of further electrons into the dot (b). 
Together with a spin-up electron from the reservoir a 
virtual spin singlet (c) is formed at energy difference 
AE = e + E12 where e is the detuning from the reser- 
voirs' Fermi energy ep- Finally the QD returns to the 
singly-charged state with a spin- up electron (d). 
E12 and Eqi are given by 



E12 — E2 — El — e 



Eoi = El — Eq — e 



(Vb - Vg) 

X 

jVg - Va) 

A 



(16) 



with Ei the energy of the QD charged with i electrons 
and A a constant describing the geometric lever arm of 
the heterostructure. 

Using (jl3p one obtains for the cotunneling rate in 
second-order [TtI. [23I 



'^cotunncl — ^rtunnol 



1 



X 



(17) 



A 



_i_ fcp 

2 tunnel 



/(e)[l-/(£)]d£ 



The integral is the sum over all second-order transi- 
tions with different detunings e from the Fermi energy 
according to Fig. [Hb-d). In addition the term with 
e{Vg — Va)/A = Eqi describes the related process where 
the virtual state is an empty QD. /(e) is the Fermi func- 
tion f{e) = 1/(1 + exp(£/fcT)). Expression dTT]) is valid 
under the condition Ez,e ^ kT i.e. for low magnetic 
fields. To obtain the exact expression for all magnetic 
fields the Fermi function terms in the integral have to be 
modified [53|]. 



The imaginary part of the denominator introduces a fi- 
nite lifetime to the electronic states limited by the tunnel- 
ing rate Ftunnei, implying that the main cause for broad- 
ening of the spin ground states is tunneling. This is rele- 
vant for elements of the integral with vanishing real part. 

In order to obtain an estimate for the cotunneling 
times in our structure we use results obtained on sam- 
ples with 25nm tunneling barrier where in certain gate 
voltage regimes tunneling rate is larger than radiative 
recombination rate i.e. Ftunnei > T, leading to broaden- 
ing in the linewidths observed in photoluminescence mea- 
surements [23^ . Then from a Wentzel-Kramers-Brillouin 
(WKB) estimation of the two different tunneling barri- 
ers together with the measured tunneling rate we esti- 
mate the tunneling rate Ftunnei to be on the order of 0.02 
to 0.1 ns~^ in our structure. We take it to be indepen- 
dent of the gate voltage within the single electron regime. 
Fig. [2] shows the calculated cotunneling rate obtained 
with expression (|17p using two different tunneling rates 
of rtunnci=0.02 ns^^ and Ftunnei=0.1 ns~^ representing 
the minimum and the maximum cotunneling rate we ex- 
pect in our experiments respectively. Co-tunneling rate is 
characterized by its very non-linear voltage dependence. 
When close to the crossover points Va and Vb it exhibits 
an ultra-steep slope; in contrast, the voltage-dependence 
is weak in the plateau center Vc- 



C. Spin - Phonon interaction 

It is known that spin relaxation in higher-dimensional 
systems is mainly due to spin-orbit (SO) interaction 
in conjunction with phonons [l^, H^- Despite being 
strongly suppressed, SO interaction is still an enabling 
mechanism for phonon-assisted spin flips in QDs and a 
considerable amount of theoretical work has been done 
on this spin relaxation mechanism [2^, [2^, [23, [H, [2^ . SO 
coupling is a well-known phenomenon in atomic physics 
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FIG. 2: Expected cotunneling rate obtained using expression 
(|17|) with the parameters Ftunnoi = 0.1 ns~^ (solid curve) and 
Ttunnci ~ 0.02 ns~^ (dashed curve), Va ~ -50 mV, Vb = 
+50mV, kT = 300 fj.eV, A=5.3. 

as well as in semiconductors and is in general character- 
ized by an interaction term of type Hso = J2i j oHjU'^j 

with I the angular momentum operator and (T the spin 
operator of the electron; the sum runs over all pairs 
i, j = X, y, z. In the case of a crystal with bulk inversion 
asymmetry (BIA) such as GaAs and InAs SO coupling 
is of Dresselhaus type [131 . Similarly Rashba SO cou- 
pling results from asymmetry along the z-direction (SIA, 
structural inversion asymmetry) [3l|. For a 2DEG the 
spin-orbit coupling can then be written as 



the piezoelectric coupling. The dependence valid for 
Ez^e. > becomes replaced by B^-kT when Ez,e < 
due to the Boltzmann factor in ([32|) , [2^, [l^ . 

In addition to this first mechanism there are ways of 
direct spin-phonon coupling which turn out to be orders 
of magnitude weaker than the admixture mechanism de- 
scribed above [1^. In [2^ the spin-flip rates due to dif- 
ferent other phonon related mechanisms are estimated 
which all depend on electronic Zeeman splitting i.e. mag- 
netic field as ~ B^. 

Beyond spin relaxation by emission or absorption of a 
single phonon, two-phonon processes in conjunction with 
SO interaction are predicted to dominate at small mag- 
netic fields. In that case a phonon with wavevector p is 
scattered into a phonon with wavevector q with energy 
conservation h/Zs\p — Ez,e- These two-phonon rates 
have characteristically strong temperature dep endence, 
estimated to be in the range T^ to T" [H, [27| . 

Alternatively it has been proposed that phonons to- 
gether with the hyperfine-induced mixing of the Zeeman 
s-levels lead to relaxation of the QD spin. As already 
mentioned in the hyperfine subsection this mechanism is 
inefficient and the resulting rate is predicted to depend 
on the external magnetic field as ~ B^ [l^ : according to 
the calculations presented in this reference the rate will 
be less than k ~ ls~^ at a magnetic field of 1 Tesla when 
considering the larger quantization energy in our QDs. 

III. GROUND-STATE OPTICAL TRANSITIONS 
OF THE SINGLY-CHARGED DOT 



Hso = P{~Px<Jx + Py<yy) + Ol{Px<Jy ~ Py^Tx) (18) 

where (3 refiects the strength of the Dresselhaus SO cou- 
pling and a the Rashba SO interaction. In ^2^ the dif- 
ference between the effect of Dresselhaus and Rashba SO 
coupling on QD spin relaxation are discussed. 

SO interaction leads to spin-orbital admixed states 
which can weakly couple to phonons leading to an ef- 
fective spin-phonon reservoir coupling of type 

^ph.off = ^^9v^-^on{hle^Qjj^^eQD,i + c.c.) (19) 
q 

Here — Cs\q\ with Cs the speed of sound. The effective 
spin-phonon coupling gphon depends on the strength of 
SO interaction, the electron-phonon coupling strength as 
well as the phonon density of states at the Zeeman energy. 

The resulting spin relaxation rate is a function of mag- 
netic field and is given by 

_ (5eAiB5)5 

where Hloq is the quantization energy for electrons and 
Ap a dimensionless constant describing the strength of 



A. Four-level model 

A singly-charged QD is described as a four-level system 
with two ground states and two excited states, coupled 
by two vertical optical transitions, as shown in FiglJlJa). 
The ground state | |) (| [)) with angular momentum 
projection = +1/2 (m^ — —1/2) is coupled to an 
excited state (trion state) formed out of two electrons 
in a singlet and a heavy hole | till) (| Ti-D-)) with spin 
projection m 2 = +3/2 (m^ = —3/2), according to optical 
selection rules by cr+ (c ) polarized optical transitions. 

The states are defined as 

I T> = f^U.+i/2lo) (21) 

I ^) = /^QD,+3/2|0> 
I Tifl") ~ ^QZ5,-l/2^QD, + l/2'''Q_D,+3/2l^) 

where Cq^, ^ {h^QD is the operator that creates an elec- 
tron (hole) in the QD with spin a along the z-axis and 
|0) is the vacuum (empty dot) state. 

All four states undergo different Zeeman shifts when 
an external DC magnetic field along the z-axis is applied. 
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FIG. 3: (a) Four-level system describing the singly-charged 
QD in magnetic field along the growth direction (z-axis). The 
electronic ground states with Zeeman splitting hcOz are ver- 
tically coupled by circularly polarized optical transitions to 
excitonic (trion) states. These consist of a heavy hole and two 
electrons forming a singlet. The fluctuations of the hyperfine 
field lead to a slowly varying coherent coupling Qh = ^H{t) 
of the spin ground states. Incoherent spin-flip processes due 
to cotunneling and phonon-SO coupling are taken into ac- 
count by relaxation rate k. A laser is introduced at Rabi 
frequency Qr and detuning Auj from the trion transition, 
(b) Transformed system after elimination of coherent coupling 
Q,H', this system is physically equivalent to that shown in (a). 
A weak hyperfine-induced diagonal transition appears at rate 
7 oc njj/B"^. The laser is now detuned on the weak 7 transi- 
tion with a reduced Rabi frequency Qr^2- 

leading to Zeeman splitting of the optical transitions. A 
cr"*" polarized laser field is introduced at Rabi frequency 
CIr and detuning Auj = ujq ~ lul with ujq the frequency 
of the trion transition and lul the laser frequency. If 
only a (7+ polarized laser field is present, the trion state 
with rriz = —3/2, i.e. (| ti-IJ-)) is inactive since the cou- 
pling strength is reduced by a factor exceeding 10'^ at 
magnetic fields larger than 60mTesla, due to a combina- 
tion of selection rules a nd, in the presence of a magnetic 
field, optical detuning [5J|. As we shall discuss shortly, 
the weak spontaneous emission to the other spin ground 
state cannot be neglected due to its long lifetime. 

Thus, the system reduces to three levels; its quantum 
dynamics is fully described by the corresponding opti- 
cal Bloch equations. These are obtained from a density 
matrix approach. 

The system Hamiltonian reads 

^ -^Zeeman ^ -^int.rad -^spin-rcscrvoir (^^) 

Then by tracing over the reservoir we obtain the master 
equation for the system (reduced) density operator p 



In the following we will discuss the different ingredients 
of this master equation. For the explicit optical Bloch 
equations we refer the reader to Appendix |B] With a 
magnetic field along the z-axis Bcxt = Bz = {0,0, B^), 
the total magnetic field at the QD is 

B = Bz+Bn (24) 

where the nuclear magnetic field (second term) is only 
seen by the electron spin, but not the hole spin. The 
Zeeman Hamiltonian then reads using (|10p 

^^Zeeman = Hz,e + HzM (25) 

Hz,e — geP-sB ■ rr 

Hz,h = ghfJ'SBz ■ (Jz 

with flnit) as defined in (fT^ . In addition, 

huj,^g,HBiBz+BNAt)) (26) 

The interaction with the radiation field in semi-classical 
form is 

Fint.rad = ftf^fl {e'^"*e^QD,-l/2^hD,+3/2 + ^'C-) (27) 

Here, the laser detuning is Aw — loq — ujl and the Rabi 
frequency is SIt?. For simplicity we use the notation 

|1) = I i> (28) 

|2) = I T) 

|3) = I Ti^) 

Hq can be written as 

/ UJ, fiHit) \ 
Ho = h\ nnit) nR (29) 

y loq-lol) 

We note here that a separation of timescales, i.e. VlH{t) ~ 
VLh can be done, since the time evolution of ^H{i) is 
much slower than all the timescales over which the system 
reaches steady-state. 

-ffspin-rcscrvoir rcfcrs to the interaction of the QD spin 
with the thermal reservoirs of electron spins and phonons. 



-^P — 7^ [Hq, p] + irelaxation (23) 

The term Hq = i/zccman + ^^int.rad dcscribcs the uni- 
tary dynamics and ^relaxation results from the interac- 
tions with reservoirs. 



-^spin-rcscrvoir — -^charge -^ph,cfF (^0) 

In our three-level model, each spin-reservoir coupling is 
treated as an incoherent relaxation rate coupling states 
|1) <-> 1 2) bidirectionally. The total rate of spin re- 
laxation, identical with the inverse of the spin Ti time 
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(k ^ = Ti), is the sum of all contributions k^. We note 
that all Ki depend on external magnetic and/or electric 
field. 

K = ^AC,; =/(yg,Sext) (31) 
— ^cotunncl(^^ J -^cxt) '^phonon(-^cxt) '^cxp 

Here, we have included a term that describes an ex- 
perimentally induced spin relaxation rate Kcxp- Unless 
specified otherwise, this relaxation is absent, but can be 
invoked by large-amplitude gate voltage modulation in 
electron cycling experiments as discussed in section ITVl 

After adding the relaxation terms due to the coupling 
to the thermal bath of radiation field modes (spontaneous 
emission terms) at rate F, the relaxation terms in the 
Lindblad form are [s^l 

^relaxation = J (2(723/0<732 ^ CTSSP — /OCTSS ) (32) 
+ f fi(2(Ti2/5(721 - 0-22/5 - /0(T22) 

-I- f (fi -I- l)(2(T2i/5o-i2 - 0-11/5- /5(7ii) 

Here, dab — is the projection operator, F is the 

spontaneous radiative decay rate of the optical tran- 
sition and K the total spin relaxation rate. At tem- 
peratures smaller or comparable to the electronic Zee- 
man splitting kT < Ez,e, a Boltzmann factor n = 
l/{exp{geiJ-BB /kT) — 1) needs to be taken into account 
which leads to thermalization of the electron spin, i.e. 
in the absence of light pw/ P22 = G^pi—Ez,e/kT), where 
Ez^e is the electronic Zeeman energy. In the case of ex- 
change coupling, the n terms cannot be regarded as an 
occupancy; it can however be shown that a similar fac- 
tor appears in the co-tunnelin g ra te (|17p when Zeeman 
splitting is taken into account 55]. 

We note, that coupling of the trion states due to hole 
spin relaxation has been neglected here. This issue is 
discussed in section ITV Gl 

B. Dressed states and rate equation description of 
spin pumping 

In order to gain a better understanding of our 3-level 
system, in the following we will transform the system 
into another basis, that gives an intuitive picture and we 
will see that the 3-level system eventually is a A-system. 
With some approximations, this allows us to capture the 
main features of the spin dynamics in the form of rate 
equations. 

The new basis is: 

|I> - |T) = cos 01 1) - sin 01 2) (33) 
It) = |2) = sin?i|l) +COS0I2) 

\m) = |3) = |3) 

with (f> = Qh/'-^z- In order to simplify the calculations 
we put the additional constraint of <C 1 and only take 



into account first-order terms in (f>. For details we refer 
to Appendix [Cl 

The transformed Hamiltonian is 

nR^2 (34) 

The off-diagonal terms due to Qh have been eliminated 
and it is obvious from the Hamiltonian ([34|) that both 
ground states couple to the excited state via an optical 
transition. Also the spontaneous emission terms become 
modified into a strong and a weak channel, marked by 
spontaneous emission rates F and 7. For details refer to 
Appendix [Cl 

The result of the transformation is shown in Figl3l(b): 
A single laser that interacted with the trion transition 
is now represented by two laser fields coupling states 1 1) 
and 1 3) (|2) and |3), respectively). Effectively, the system 
can be decomposed into two two-level systems with its 
own spontaneous emission rates, Rabi frequencies and 
effective laser detunings. Those are for the |1) ^ |3) 
subsystem: 

^ 7 = (35) 

^R,l = i'^R 
AtJi — Aw -I- LOz 

And for the |2) ^ |3) subsystem: 

F = F (36) 

^R,2 — ^R 
AlU2 ~ Auj 

It is clear that, in this first order approximation, the 
couplings of the dressed |2) ^ |3) system are the same 
as for the bare |2) ^ |3) transition. 

In the following we will discuss the properties of the 
transformed three-level system with a resonant laser on 
the 1 2) ^ 1 3) subsystem, i.e. under the condition Aw = 
0. 

For that, we start out with the system being in state 
1 2) and we ignore state |1). After an intermediate time 
to given by ^ to ^ the laser field induces a 

steady-state occupation of the excited state ps^ito). Fur- 
ther, for times much longer than (p33(to)7)^^ the system 
can also be found in state |1). 

The net effect of this spin-flip Raman process is a trans- 
fer of occupation from state |2) to state |1). We will refer 
to this process as OSP, due to its similarity to experi- 
ments performed with atoms [33| . Further, we note that 
a scheme that uses OSP for spin state preparation had 
been proposed in [s^l, and state preparation has been 
experimentally demonstrated in Q. 

In order to obtain the OSP rate that transfers the sys- 
tem from state |2) to state |1) under the presence of a 
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resonant laser {Acu = 0), we assume the system is in 
state 1 2). As already mentioned, 7 represents a weak es- 
cape channel only, and we can work with the steady state 
occupation for the trion state |3), which is [ssj 



^2 



Y2 



(37) 



The OSP rate from state |2) to |1) then is the trion steady 
state occupation times the spontaneous emission rate into 
the weakly allowed channel: 



R2 



00) ■ 7 



(38) 



The |1) <-> 1 3) sub-system has a spontaneous emission 
rate of 7, however the relaxation of its excited state is 
governed by a strong escape channel which is determined 
by the F rate. The rate of transfer from |1) to |2) is ^ 



Ri- 



02 r 

4tj? + f 2 



n2 r 

4w? 



(39) 



Using ([55]) to ([M]) . the time-averaged ratio of ground 
state occupations is obtained under the condition that 



P22 
Pii 



R2^1 



2n% 



4co? 



(40) 



The approximation on the right hand side is valid in the 
limit of a weak incident beam (fi^j <C F) and a Zeeman 
splitting largely exceeding the trion decay rate {uj^ ^ F). 
Under these conditions, from ([57)1 follows that p33(oo) ^ 
1, such that together with tr(p) = 1 we can safely assume 



Pii + /O22 « 1 



and we obtain 



P22{t = 00) 



1 



(41) 



(42) 



r2 



In the case i?i^2 ^ k, i?2^i the ratio of ground state 
occupations is 



P22 



Pll R2^1 + K 

and together with (|¥T|) we obtain 

1 

P22(t = 00) 



C 



2 + C 

7 nl 



(43) 



(44) 



K F2 -f 2n% 



Hence in the case of fixed laser intensity i.e. constant flj^, 
the spin-state occupations are determined by the ratio of 
OSP rate versus spin relaxation rate. 



C. Hole mixing 

Valence-band mixing, as described by the Luttinger 
Hamiltonian (36j . is a well-known feature in quantum 
wells. Similarly, it is expected to play a role in quantum 
dots. With valence-band mixing, a heavy hole acquires a 
small contribution of light holes and vice versa such that 
the effective hole state as it was defined in pij) has the 
form 



(4 



QD,+3/2 



(45) 



with |£±| <C 1. Pscudopotcntial calculations for self- 
assembled InAs QDs yield admixtures on the order of 
a few percent [37]. As it has been pointed out in Q 
valence-band mixing would have a major impact upon 
the effective optical selection rules by introducing a diag- 
onal relaxation channel between states |3) and |1) due to 
the admixed light hole component of state |3). Two cases 
have to be distinguished: First, the mixing contribution 
associated with e+ further leads to an effective coher- 
ent laser coupling in addition to the coupling induced by 
hyperfine interaction (cx ^j^^^fl) at a detuning Aw -I- tOz 
as shown in Fig. [3jb). Second, the e_ part essentially 
only appears as a relaxation channel without coherent 
laser coupling, as the dipole moment of this linearly po- 
larized transition lies along the propagation axis of the 
laser beam and therefore cannot be excited. Hence the 
following diagonal relaxation terms are added to equation 
(El 



* _ 7hm - 

-t^rolax.hm — ['^0'13 P<^31 



0-33/0 - P^33) (46) 



with7h,n= |e|2r= (|e+|2 + |e_|2)r. 

This diagonal rate leads to OSP in a way similar to the 
7 channel enabled by hyperfine interaction. The main 
difference between hyperfine-induced OSP and valence- 
band mixing induced OSP is that the first one is magnetic 
field dependent as discussed previously, and the latter is 
not: the valence-band mixing strength e is expected to 
be independent of magnetic field, as long as the Zeeman 
splitting is much smaller than the heavy-light hole split- 
ting (A/i/ >10meV), which is true for all realistic experi- 
mental magnetic fields. Since the hyperfine-induced OSP 
rate drops with magnetic field 7 cx B~^, hole- mixing in- 
duced OSP should dominate at high fields. From our 
measurements at high magnetic fields (FiglHJb)) we ex- 
tract a 7^^^ of 2 ± 0.8/xs which yields a hole mixing 
strength of \e\ ^ 2.2%. This |e| value is indeed much 
smaller than 1, but we also expect the exact value to 
vary from one QD to another. 

We note that a slightly tilted external magnetic field 
would yield identical dynamics in the absence of any hole 
mixing since it would lead to mixing of electronic states 
induced by the in-plane component of the applied field. 
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These two fundamentally different mechanisms are ex- 
perimentally indistinguishable for a fixed magnetic field 
orientation. Therefore we repeated our experiments as 
a function of sample tilt under a magnetic field. For a 
±1.5° coverage of tilt in all directions our measurements 
yielded no observable change in the measured quantity 
"fhm- Hence, we can safely state that the inherent hole 
mixing in our QDs indeed dominates over small-angle 
tilt-induced mixing of electronic spin states. 

IV. SINGLE DOT ABSORPTION 
SPECTROSCOPY WITH RESONANT LASER 

A. Experimental method 



the optical transition, u the laser frequency and v the 
scattering volume of the QD. The real part of x leads to 
a dispersive QD response to the laser detuning whereas 
the imaginary part produces an absorptive Lorentzian 
response as can be read from (|47p . 

The signal detected in a DT experiment arises from 
interference of the forward scattered field together with 
the excitation field jiQj. This can also be seen from the 
optical theorem p[3 | which relates the absorption cross- 
section of the dipole to the forward-scattering amplitude. 
When the QD is exactly at the focus, the imaginary part 
of the scattered field is in phase with the excitation laser. 

After collecting all factors describing spatial mismatch 
between excitation field and QD scattering cross-section 
we define relative absorption as 



All data shown in this work has been obtained u sing 
differential transmission (DT) technique [H, [H, HO, |4l[ 
|42 . 43] . A narrow-band laser is scanned over the QD res- 
onance while the transmitted light intensity is measured 
by photodetectors (for details refer to Appendix [X)) . A 
QD transition then results in a transmission dip i.e. ab- 
sorption due to resonant Rayleigh scattering. In ^ the 
DT technique is used to investigate fine structure of the 
X° which exhibits a single two-level system in the case 
that the scanning laser is linearly polarized along one of 
the QD axes and only one transition of the AT^ doublet 
is addressed. 

In order to link the experimentally observable absorp- 
tion i.e. intensity drop on the detectors to the 3-level 
system of the singly- charged QD, we use the effective A- 
system picture as described in the previous section. With 
a resonant laser and large external magnetic field along 
the z-axis i.e. Aw = 0, F ^ 7, the 2^3 sub- 

system with strong spontaneous emission rate F acts as 
the main scattering source. In the B ^ 00 limit only the 
strong transition contributes to light scattering. 

The QD response in this type of DT experiments is dis- 
cussed in the above mentioned references [iO, |4l| . In the 
following we therefore only sketch the link between our 
3-level system and the intensity of the light transmitted 
through the sample. 

In general, for a given two-level system with ground 
state \g), excited state |e) and spontaneous emission rate 
F the linear susceptibility x — Xqd defined by P = 
EqExE, with P the polarization describing the QD dipole, 
E the electric field, Eq the dielectric constant of vacuum 
and e the dielectric constant of medium (GaAs) . x can be 
written in terms of the steady-state off-diagonal density 
matrix element peg{oo) 



x(^) 



2Trc^v 



f (c^o-cj + iF/2) 



(47) 



Wo' (wo - w)^ + + "1/2 

2nc^v F , , 
= 

with c the speed of light in vacuum, ujo the frequency of 



e(Aw) = 1 - 



T{Au 

Tnft 



= s.Im(-— pe5(oo)) (48) 



where T^Au) refers to the transmitted intensity as a 
function of laser detuning and Toff to the transmitted 
intensity far off resonance in the limit (Aw — > 00). On 
resonance in the weak excitation regime i.e. fi^ <C F 
the term Im( jy^peg(cx))) reaches 1 and s is a scaling fac- 
tor characterizing the maximum theoretical absorption 
contrast that is given by 



(49) 



valid for a weak focusing geometry. Here ao — 
(3/27r)(A/n)^ is the scattering cross-section of the two- 
level system in the weak excitation limit and Al is the 
laser spot area. The factor S'exp accounts for reduction 
of signal due to our lock-in detection scheme and exper- 
imental imperfections. Further, with Aw — 



Peg{oo) = -i— — Pee(oo) 



fir 



(50) 



and equation (|^ reads 



9(Aw = 0) = 1 - ^^^^ Q) ^,.l!-p^^(oo) (51) 

J off "fl; 

We can now apply this scheme to our double 2-level sys- 
tem and we will assume in the following that the laser 
is^ on resonance with the strong trion transition i.e. the 
1 2) <-> 1 3) subsystem. Clearly, ^33 = pee and P23 = Pge- 
When llIz S> the weak transition does not contribute 
to the absorption signal. Using ((5T|) we can infer the 
value of the spin-up state occupation P22 (00) from our ab- 
sorption measurements when varying parameters such as 
magnetic field and gate voltage but keeping laser power 
constant i.e. constant il/j. 

We still need a calibration point, i.e. an experimental 
value of 6(0) for a known (022(00). In the absence of 
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an external magnetic field the spin ground states can be 
considered to be fully mixed due to the in-plane part of 
the Overhauser field leading to 7 ~ F and a^branching 
ratio of 77 = 1. As a consequence the |1)-|3) and the 
|2)-|3) transitions equally contribute to light scattering 
and fast bidirectional OSP takes place, leading to a fully 
randomized spin i.e. 



probe laser detuning has been readjusted in the graphs 
to compensate the Zeeman splitting. 

Absorption drops by nearly two orders of magnitude 
over the plotted range of — to 300 mT. With a 
resonant laser in the weak excitation limit and Zeeman 
splitting much larger than the trion transition linewidth 
i.e. Au! — 0, 3> r ^ 7 equation yielded 



pii{t = oo,B = 0) = p22it = oo,B^ 0) 



(52) 



Given the relations ((5T|) and ((52|) and the fact that the 
absorption measurement constitutes a relative measure- 
ment of the steady-state occupation ^22(00) and expres- 
sion (jSip can be rewritten by replacing the s-factor 



e(Aw = 0) 



^ _ TjAu; = 0) ^ ^, 

Toff ft 



r2 ^ 

7^P22(oo) 



(53) 



The s'-factor can be experimentally determined using 

(El). 



B. Optical spin pumping 



(a) 



(b) 
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FIG. 4; (Color online) (a) Absorption maxima in the plateau 
center plotted as a function of magnetic field B^. A drop 
occurs with increasing due to OSP which, at low mag- 
netic fields, dominates over co-tunneling and phonon inter- 
action. The grey line is a numerical simulation using Q = 
0.6 r, _Bnuc = 15mT, r~^=0.8ns, additional diagonal relaxation 
Thm ~ 2/^s, K~^ = 10ms. The inset shows the corresponding 
raw laser scans from T (top) to 300 mT (bottom). The peaks 
have been shifted laterally to eliminate Zeeman splitting, (b) 
Optically induced spin pumping rates -R2->i, transferring the 
system into the dark state, and Ri^2 the back-pumping rate. 

FigHUa) shows absorption on resonance on the blue 
Zeeman transition as a function of magnetic field nor- 
malized to on-resonance absorption at Tesla i.e. 
0(Soxt)/0(5oxt = 0). The gate voltage was kept in the 
plateau center i.e. in a regime where Kcotunnoi is minimal. 
The inset shows the corresponding raw laser scans for 
T (top) to 300 mT (bottom). The zero positions of the 



P22(t = 00) 



1 + ^ 



all provided that the spin relaxation rate k <C 
i?2— »i,^^i^2 which we can safely assume for low mag- 
netic fields [H, and strongly suppressed exchange 
coupling in the gate voltage plateau center. Consistent 
with (|42| . drop of absorption follows a B^^ cx law 
indicated by the dashed line. For fields less than lOOmT 
(see Fig. IH^a)) the approximations included in (|^^ do 
not hold any more and F 7. Without any approxima- 
tion the steady-state solutions of the optical Bloch equa- 
tions are evaluated (solid line) numerically; they are in 
excellent agreement with our data at all magnetic fields. 
The Rabi frequency f2 in units of F for a given incident 
laser power can be independently determined by satura- 
tion spectroscopy and power broadening measurements. 
The radiative lifetime F~^= 0.8ns used in our simulation 
is based on a measurement in as-grown dots [45| . 



C. Electron cycling 

Given that the OSP rates Ri^2, -^1^2 and the spin re- 
laxation rate k are unknown the experimental data shown 
in Fig. HJa) do not reveal direct quantitative information 
about 7. However, the branching ratio can be extracted 
using an rms-coherent coupling {n'j^{t)) given in (|Dip 



7 _ {nj,{t)) Bl 



F-h7 



25? 



(54) 



?7 is equivalent to the probability that the system decays 
via the 7-channel when excited into a trion state. 

To determine rj we applied a large square- wave modu- 
lation (Amplitude 80mV peak-peak) at different frequen- 
cies to the gate which in every cycle first loaded another 
electron of opposite spin into the QD forming a singlet 
together with the QD electron. Then one of the electrons 
was forced to leave, and as the tunneling probability for 
each of the two electrons is equal the remaining QD spin 
was fully randomized. The advantage of this technique 
which we will refer to as electron cycling leads to en- 
forced spin relaxation at a known rate Koxp (also refer 
to equation (|3ip ). In the case k « Koxp ^ Ri^2 i-e. 
enforced spin relaxation rate exceeds the optical back- 
pumping rate i?2->i and 7 can be determined by a fit 
using equations (|43l44p . 
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FIG. 5: (Color online) (a) Electron recycling measurements. 
Peak absorption in the plateau center normalized to peak ab- 
sorption in the co-tunneling regime is plotted as a function 
of laser power at a constant _Bz=300mT and three different 
K ~ Kcxp=54 kHz (upper, green points), 19kHz (middle, red 
points) and 3kHz (lower, blue points). The fits indicated by 
the solid grey lines have been obtained using equation (|44|) 
yielding 7^^= 0.63 fis. (b) A check experiment: Two laser 
scans at K = Kcxp=54 kHz with in-and-out of plateau modula- 
tion (showing peak) and with in-plateau modulation, demon- 
strating that indeed controlled spin relaxation is realized. The 
noise level is indicated by the dashed blue line, (c) Intensity 
dependence of relative absorption at B^^OT (red circles) and 
-Bz=100 mT (blue squares). There is essentially no depen- 
dence on laser power confirming the theoretical model which 
gives p2|l . 



Fig. EJa) shows absorption normalized to on-resonance 
absorption in the co-tunnehng regime; the data was ob- 
tained with electron cycling for different modulation fre- 
quencies and laser powers at a fixed external magnetic 
field Bz = 300 mT. The upper, green points correspond 
to K=54.3kHz, the middle, red point to k =19.3 kHz, and 
the lower, blue points to k=3.3 kHz. 

Using ([5^ . (|5ip and (I44p the absorption ratio shown 
in the Figure can be written as 



^inplateau /^22,inplatcau 



^cotunnel 



P22,cotu 



2 + C 



(55) 



The grey lines are fits using this expression. Best match 
with the data can be obtained with a total OSP rate 
7tot(300mT)=7(300mT)-|-7,,,„= 1.6^is^^ [57], where the 
two contributions stem from nuclear spins and hole mix- 
ing, respectively. The hole mixing-induced contribution 
can be independently determined from high-magnetic 
field measurements to be "fhm = 2^s. Using the branch- 
ing ratio ([51)1 with F^^ = 0.8ns we then solve for the rms- 
nuclear magnetic field and obtain -Bnuc = 15mT which is 
in good agreement with 



Fig. mb) shows a measurement that demonstrates 
the difference between electron cycling (large amplitude 
modulation) and in-plateau (small amplitude) modula- 
tion: For in-plateau modulation we do not observe ab- 
sorption (the noise level is marked by the horizontal 
dashed line). In contrast, when large amplitude mod- 
ulation is applied absorption is partially recovered due 
to forced spin relaxation at a controlled rate Kexp, as 
shown by the red peak in the figure. Figl^c) was ob- 
tained in the plateau center without electron cycling tech- 
nique, showing on-resonance absorption as a function of 
incident laser power. At 100 mT (blue squares) relative 
absorption is one order of magnitude weaker than at 
mT (red circles) due to OSP. In both cases absorption ex- 
hibits weak dependence on laser power. This dependence 
arises from being in the vicinity of laser power required 
for saturation, hence a deviation from the assumptions 
of 



D. Peak shift in the plateau center 
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FIG. 6: Gray-scale plot of absorption strength as a function 
of laser detuning and gate voltage. External magnetic field of 
-Bz=150mT is applied. The linear gate- voltage dependence is 
due to the quantum-confined Stark effect, (a) and (c) show 
the absorption in the vicinity of the blue optical transition, (b) 
and (d) show the red transition, (a) and (b) are experimen- 
tal data showing a spectral shift of the weak absorption peak 
in the plateau center compared to the strong co-tunneling 
regimes at the plateau edges. The shift is directed to the 
lower energies for the blue transition, and towards higher en- 
ergies for the red transition. This feature can be reproduced 
in a numerical simulation including a randomly fluctuating 
Overhauser field as shown in (c) and (d) using the parame- 
ters n = 0.6 r, r-^=0.8ns, Snuc=15mT, 7-^=2/is. 

Fig. [S] shows laser scans obtained at i32=150mT on 
the blue (a) and red (b) Zeeman transition throughout 
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the whole single-electron plateau. Absorption strength is 
grayscale-coded. The line tilt is due to the quantum con- 
fined Stark effect and pixelization is due to experimen- 
tally limited voltage resolution. At gate voltages 395 mV 
and 480 mV the two co-tunneling regimes show strong ab- 
sorption when spin relaxation is fast due to charge reser- 
voir coupling (for details refer to IIIB|) . In the plateau 
center hyperfine interaction dominates and leads to spin 
pumping and drop of absorption as already discussed. 
Here, we further observe a shift of the spectral position 
of the absorption peak in the pleateau center as com- 
pared to the co-tunneling regime. This shift is directed 
to the red (blue) for the blue (red) Zeeman transition. 
This resembles effects one might expect for dynamical 
nuclear spin polarization (DNSP) [io, [13]; these effects 
can however be excluded 58 ]. Our numerical simulation 
is able to reproduce this behaviour without taking into 
account DNSP as shown in Fig. [H] (c) and (d). The 
lineshift in the plateau center is ±0.9GHz for both red 
and blue transition which is close to the electronic Zee- 
man splitting at 150 mT, Ez^e ~ 1-3 GHz. At the co- 
tunneling edges K is large and maximum of absorption 
is observed when the laser is exactly on resonance with 
the transition. When k, is small as it is the case in the 
plateau center absorption of a strictly resonant laser is 
suppressed due to spin pumping depending on the exter- 
nal magnetic field. When the laser frequency is moved 
towards the center between the strong T and the weak 7 
transitions i.e. the spectral detuning with respect to the 
7 transition is reduced, the back-pumping at rate Ri—,2 
becomes more efficient and maximum of absorption will 
be reached for a spectral detuning that fulfils the condi- 
tion Ri^2 — R2—>i- As a consequence both transitions 
contribute to absorption which leads to a shift of the ab- 
sorption maximum towards the weak 7 transition i.e. a 
blueshift when the red line is observed and vice versa. 



E. Peak broadening at plateau edges 

Fig. [7{a) shows example laser scans for different mag- 
netic fields Bz ranging from to 1 Tesla, obtained in the 
co-tunneling regime where spin relaxation is fast. The 
scans have been laterally shifted in order to eliminate 
the Zeeman shift. In (b) the measured linewidths are 
plotted as a function of magnetic field (red circles) along 
with a calculated curve (solid line). A broadening to 
almost double the zero-field linewidth appears at mag- 
netic fields between 60 mT and 80 mT; at higher fields 
linewidth becomes as narrow as in the case = 0. 

The physical reason for the observed broadening is 
very similar to that described in the previous paragraph 
IIVDI Both 7 and F transitions contribute in a non- 
negligible way to absorption and maximum is observed 
when Ri-,2 = R2^i condition is fulfilled. Consistently 
the linewidth increases as much as the electronic Zee- 
man splitting initially but drops at magnetic fields where 
T{B) 3> ^{B) and a single transition is established. 




Laser Detuning (GHz) Magnetic Field (ml) 

FIG. 7: (Color online) (a) Trion transition laser scans for 
five different magnetic fields. The gate voltage was in the co- 
tunneling regime (see section fll Bp . (b) Measured linewidth 
obtained from laser scans as a function of magnetic field indi- 
cated by the red circles. A broadening occurs at Bz ~75 mT; 
at larger fields the linewidth almost recovers back its original 
value of 450 MHz at Tesla. The solid lines are obtained via 
numerical simulation with parameters fl — O.GF, i3nuc=15mT, 
r-^=0.8ns, 7-^ = 2.0AIS, k-^=2.5ms. 

In contrast to section IIVDI k is large here due to co- 
tunneling and annihilates the absorption drop caused by 
spin pumping hence making the transition visible at all 
magnetic fields. The solid lines are calculated curves for 
a range of e values using a randomly fluctuating Over- 
hauser field with -Bnuc = 15mT, well reproducing this fea- 
ture. We note that in order to put as many constraints 
as possible on the choice of simulation parameters we 
have used the maximum co-tunneling-induced spin re- 
laxation rate k = Kcotunnoi = OAfis~^ as obtained from 
the data shown in FigEfb). The effect of e on the simu- 
lation is clearly negligible, advocating that the dominant 
OSP mechanism is hyperfine interaction. 



F. Coupling to electron spin reservoir 

We have performed laser scans as a function of gate 
voltage througout the whole single-electron plateau as 
defined in section III B I The measured on-resonance ab- 
sorption signal for each laser scan is plotted in Fig. [S^a). 
The data has been obtained at an external magnetic field 
of -82=300 mT taking coarse voltage steps. 

At gate voltages lower than 540 mV and higher 
than 625 mV as marked by the shaded regions absorp- 
tion drops below noise level indicated by the horizontal 
dashed line. At these voltages the QD either becomes 
empty (left of point A) or doubly charged (right of point 
B). Absorption then vanishes since in those cases the QD 
is not described by the trion level system any more; the 
optical transitions for these gate voltages are not ob- 
served within our scanning window of 30 GHz around 
the trion transitions. The unshaded part indicates the 
region where the QD contains a single electron and as it 
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ation mechanism. The gray sohd Une is a best-fit numer- 
ical simulation using expression (|17p as spin relaxation 
rate showing good accordance with the data. The co- 
tunneling rate at the peak as determined from the fit is 
'^max~2.5/is. The noise level is indicated by the dashed 
line; it deviates from the one shown in (a) due to different 
experimental settings such as lock-in time constants and 
filter slopes. Again on the left side of the peak the QD 
is empty, yielding vanishing absorption below the noise 
level. The gradual decrease of absorption is due to finite 
temperature. On the right side the spin pumping regime 
is located; here some weak absorption remains according 
to the occupation of the observed spin state, revealing 
the strength of spin pumping. 



FIG. 8: (Color online) (a) Example voltage coarse scan across 
the entire single-electron plateau at _B2=300mT. Per voltage 
step a laser scan is performed and the observed absorption 
maximum is plotted as a function of gate voltage. Due to 
strong hyperfine-induced OSP and weak co-tunneling rate ab- 
sorption in the plateau center is suppressed; however when 
approaching the single-electron plateau edges absorption is 
recovered due to highly nonlinear dependence of co-tunneling 
on gate voltage leading to fast spin flips. Outside the volt- 
age plateau i.e. left of point A or right of point B absorption 
is suppressed because the QD then becomes either empty or 
doubly-charged which shifts the optical transition energies out 
of our spectral observation window of 30GHz. The solid line 
is a guide to the eye. (b) shows a voltage fine scan of the 
left co-tunneling regime obtained on another QD than in (a) . 
The solid line is a numerical calculation using F"^ — 0.8ns, 
Q=0.6 F, Bnuc =15mT, r-J„„^i=20ns, 7"^ = 2fis and a Zee- 
man splitting of Ez,e = 10/xeV. The voltage-FWHM of the 
co-tunneling peak is lOmV. 



G. Coupling to phonon reservoir 




Gate Voltage (mV) Magnetic Field (T) 



has been mentioned before the co-tunneling rate is max- 
imum when gate voltage is at the crossover points A or 
B. Here, relaxation via co-tunneling is faster than the 
optical pumping rates Kcotunnoi ^ Ri~>2,R2~*i leading 
to thermalization of the electron spin and thus strong 
absorption. The scenario drastically changes when gate 
voltage is tuned to the center of the plateau. Here, co- 
tunneling rate Kcotunnei rcaches its minimum where our 
numerical calculation predicts a drop of as much as five 
orders of magnitude (also see Fig. ^ compared to the 
crossover points such that Kcotunnoi ^ R2—fi- Conse- 
quently the occupation of the spin states is governed by 
OSP (equation (|42)) ) rather than Boltzmann factor mean- 
ing that the spin is predominantly in the dark state and 
vanishing absorption is observed. 

The semilogarithmic plot in Fig. [D^b) shows a volt- 
age fine scan of the low voltage plateau edge around 
the A crossover point obtained at i?z=300mT. The gate 
voltage for point A is different from Fig. ^a,) as this 
data was taken on another QD. The observed absorption 
drops by half within a gate voltage detuning of ±5mV 
from the maximum position. This data demonstrates 
the enormous gate- voltage dependence of this spin relax- 



FIG. 9: (Color online) (a) Plateau laser scans for 4 different 
magnetic fields: Absorption is plotted as a function of laser 
detuning and gate voltage. Zeeman effect has been eliminated 
by vertical shift of each single plot. The linear voltage depen- 
dence is due to quantum-confined Stark effect. Plateau center 
absorption drops at intermediate magnetic fields due to OSP 
but then recovers at high fields due to fast thermalization via 
phonon-SO interaction at 9.9 Tesla. At the plateau edges the 
QD can absorb at all magnetic fields due to fast /^cotunnci . 
(b) Black data points: Magnetic-field evolution of peak ab- 
sorption in the plateau center normalized to peak absorp- 
tion in the co-tunneling regime. The lower (upper) boundary 
of the grey region is obtained from a numerical simulation 
for 7,7^ = 1.2/^s (7^^ = 2.8/^s); the red line corresponds to 
Ihrn ~ 2/^8. Bcst match with the data yields explicitly for the 
spin relaxation rate Kphonon = aoB^ with ao — 0.031 in units 
of [T"^sec"^]. Again, Bnuc = 15mT. 

Based on the theoretical estimates of section III CI we 
now seek for signatures of SO-phonon induced spin relax- 
ation at high magnetic fields. Fig[5]Ja) shows 2D plots of 
color-coded absorption strength as a function of laser de- 
tuning and gate voltage for four different magnetic fields 
obtained for the red Zeeman transition. The scans cover 
the whole single-electron plateau; excitonic Zeeman shift 
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has been eliminated by shifting the y-scale for each 2D 
graph separately. The linear dependence of the excitonic 
transition energy on gate voltage is due to the quantum 
confined stark shift. 

At Tesla absorption is clearly visible throughout the 
whole plateau due to fast spin flips with the neighboring 
nuclear spins. When a small magnetic field (B=0.1T) is 
applied, absorption in the plateau center drops because 
of hyperfinc-induccd OSP as disucsscd in the previous 
sections. Close to the plateau edges absorption still re- 
mains due to fast co-tunneling. At 0.5 Tesla increasing 
OSP leads to further drop of absorption. These absorp- 
tion characteristics in the plateau center remain the same 
up to 5 Tesla, however absorption starts to come back at 
even higher fields: when the magnetic field is raised up to 
9.9 Tesla a significant recovery of plateau-center absorp- 
tion is observed. This effect cannot be explained by OSP 
which only causes monotonous decrease of absorption nor 
by co-tunneling which is negligible in the plateau center 
and hardly shows any magnetic field dependence. Ow- 
ing to its dependence however, phonon-assisted spin 
relaxation is a good candidate to be responsible for the 
observed effects in the context of spin relaxation mecha- 
nisms. 

FiglSlJb) shows the quantitative evolution of normal- 
ized absorption with magnetic field i.e. the ratio of ab- 
sorption in the plateau center versus co-tunneling regime 
(black data points). Further, the solid red line along 
with the gray shaded region indicate calculated strength 
of absorption for = 2/zs with an uncertainty of 

±0.8/j,s; the phonon-induced spin relaxation rate was 
Kphonon = C(oB^ with the coefficient ao = 0.031 in 
units of [T^^sec"^]. Whereas Kphonon is strongly B- 
dependent, the hole mixing contribution has no B- 
field dependence within the magnetic field range consid- 
ered here. Therefore these two mechanisms have distin- 
guishable effects on FiglHlJb) and thus can be identified in- 
dependently. The good agreement with the experimental 
data strongly suggests that in this regime of electric and 
magnetic fields the dominant spin relaxation is indeed 
phonon-assisted. Further, within our uncertainty Kphonon 
matches well with the results that have been previously 
obtained on an ensemble of self-assembled InAs/GaAs 
QDs 11]. 

There are two fundamentally different mechanisms 
which employ holes to yield OSP: First, hole mixing of 
strength e leads to an admixture of the light hole states 
to the trion states as dicussed in section IIII CI Second, 
hole-spin relaxation leads to an incoherent coupling of 
the trion states contributing to OSP. In 48] hole spin re- 
laxation rate is predicted to be below lO'^/sec and mono- 
tonically increase with magnetic field which suggests that 
it is not the main mechanism responsible for OSP. We 
therefore neglect hole spin- flips, further assuming that 
there is no other efficient hole spin flip mechanism at low 
magnetic fields. In the first mechanism, OSP is indepen- 
dent of magnetic field and the strength is equal to the 
hyperfine- induced OSP rate at '--^l Tesla. At higher mag- 



netic fields the hyperfine-induced OSP rate drops with 
B~^, therefore hole mixing becomes the dominant OSP 
mechanism here. 
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FIG. 10: (Color online) Calculated absorption maxima for 
the whole single-electron plateau plotted as a function of mag- 
netic and electric field, (a) shows the simulation for the probe 
laser in the vicinity of the red Zeeman transition, (b) simi- 
lar but for the blue Zeeman transition. The borders of each 
plot show strong absorption due to interactions with nuclear 
spins (left), charge reservoir (top/bottom) and phonon reser- 
voir (right). At large magnetic fields spin polarization nearly 
reaches unity due to thermalization leading to vanishing ab- 
sorption on the red transition (a), and enhanced absorption 
on the blue transition (b). In the center of the plots absorp- 
tion and thus spin relaxation is suppressed by approximately 
five orders of magnitude. The parameters used in the simula- 
tion are; r^^=0.8ns, _Bnuc = 15mT, j^^=2fis, tunneling time 



=20ns, K, 



, ^phonon 



as given in section HV CI 



In the previous sections the three dominant mecha- 
nisms acting on the confined spin have been identified 
separately along with the two mechanisms for OSP. In 
this final part we will present results of our numerical 
simulation based on parameters that have been measured 
or estimated before. The calculations have been per- 
formed within a parameter space approximately overlap- 
ping with the full scale of our experimental tuning ability 
of the static electric and magnetic fields. For details of 
the simulation we refer to Appendix [Dl 

Fig. [ini shows calculated maximum values of absorp- 
tion for the red (a) and the blue (b) trion transition with 
a laser having the corresponding circular polarization. 
Absorption strength is color-coded in logarithmic scale 
as a function of gate voltage detuning and external mag- 
netic field. 

All of the following points have been discussed in the 
previous chapters; here we mention them briefly as a 
key to the plots: the necessary conditions for observ- 
ing strong absorption are either k ^ i?i_,2,-R2^i or 
i?i^2 ~ -^2^1- Further, at large magnetic fields B > 8 
Tesla when Ez.e ^ kT the Boltzmann factor leads to a 
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difference of the spin ground state occupations and thus 
a difference between absorption strength on the red and 
blue Zeeman transition. 

In the plot we distinguish three different regimes of 
strong absorption: 

1. Magnetic fields lower than the fluctuations of the 
hyperfine field (B ^15 mTesla). Here, fast bidirec- 
tional OSP due to hypcrfine-induced state mixing 
leads to strong absorption. 



Here, '^phonon 
^phonon 



2. High magnetic fields (>5 Tesla) 
induces fast thermalization i.e. 
i?i_»2, i?2— >!• The spin ground state occupation is 
mainly determined by the Boltzmann factor leading 
to a lowering (increase) of absorption on the higher 
(lower) energy spin state occupation (a) ((b)). 

3. Large gate voltage detunings from the plateau cen- 
ter (±40mV). Here, co-tunneling (Kcotunnei) is re- 
sponsible for fast spin relaxation and appearance of 
absorption. 

An intriguing feature that becomes apparent now is 
the blue island in the center of the color-coded spin re- 
laxation plot. It marks the regime where absorption (i.e. 
all reservoir interactions) is suppressed by five orders of 
magnitude, or in other words the localized spin becomes 
maximally isolated hence the frequently used concept of 
an artificial atom is meaningful. Within the scope of 
quantum information processing this indicates the rele- 
vant regime of operation. 



VI. SUMMARY AND CONCLUSION 

We have investigated the dominant interactions of a 
confined electron spin in a single self-assembled QD by 
optical means and demonstrated the regimes where each 
reservoir coupling becomes important. For magnetic 
fields B < ITesla, the dominant contribution to OSP 
stems from the fluctuating hyperfine field mixing the elec- 
tronic spin states and creating a weak channel for diag- 
onal relaxation in the trion four-level picture. Exchange 
and phonon-induced spin-flip processes dominate over 
hyperfine-induced spin pumping and establish a thermal 
steady-state at large gate- voltage detunings and/or large 
external magnetic fields; in the plateau center at inter- 
mediate magnetic fields the situation is reversed and spin 
pumping dominates, strongly altering the state occupa- 
tions away from thermal equilibrium values. Signatures 
of heavy-light hole mixing dominated spin cooling can be 
observed for fields > 5Tesla. 

From a quantum control perspective these results 
demonstrate that the quantum dynamics of a single con- 
fined spin can be significantly altered by externally con- 
trolled parameters such as electric and magnetic fields. 
A natural extension of this study would be the investiga- 
tion of spin decoherence in a single QD using similar opti- 
cal techniques. These measurements would require more 



advanced schemes such as electromagnetically-induced 
transparency (EIT). Further, knowledge gained on single- 
electron spin dynamics can be utilized in the resonant 
optical study of more complex systems such as coupled 
QDs or QDs with a single excess heavy-hole. 
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APPENDIX A: SAMPLE AND EXPERIMENTAL 
TECHNIQUES 
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FIG. 11: (a) Layers of the Schottky-type heterostructure. (b) 
Corresponding band structure. At gate voltage Vb the lowest 
QD conduction band level lies above the Fermi energy and 
the QD is empty, (c) At gate voltage Vi > Vb the lowest QD 
level is below the Fermi energy and therefore populated with 
one electron. A second electron cannot enter due to required 
charging energy {Coulomb blockade). 

Our InAs/GaAs quantum dots (QDs) are grown 
by molecular beam epitaxy in Stranski-Krastanow 
mode leading to lens-shaped dots of average size 
25nmx25nmx5nm; QD light emission is blue-shifted by 
partially covered islands (PCI) technique. A 35-nm GaAs 
tunneling barrier separates the QDs from a charge reser- 
voir formed by a heavily doped n-GaAs layer which forms 
the back contact. Above the QDs there is a 12 nm thick 
GaAs cap and a 50 nm-thick Alo.4Gao.6As blocking layer 
which prevents the holes from coupling to the continuum 
states within the 88-nm capping layer [22|. Bias voltage 
between the back contact and a semi-transparent 5 nm- 
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Ti-Schottky window determines the electric field in the 
structure and allows us to load a single conduction-band 
electron into the QD. 




Vg (mV) 



only apply resonant excitation of the QD single electron 
ground state transitions by using a differential transmis- 
sion (DT) technique O, EH, [H [11] . In order to obtain 
a spectrum we sweep a single-mode Ti;Sa laser over the 
QD transition and record the intensity of the transmit- 
ted light. A QD resonance is observed as a dip on top of 
the laser background. The resolution of this technique is 
only limited by the laser linewidth i.e. /S.VLaser <1 MHz. 



APPENDIX B: OPTICAL BLOCK EQUATIONS 
FOR THE THREE-LEVEL SYSTEM 

The optical Bloch equations are derived from the mas- 
ter equation (|23p . Including rotating- wave approxima- 
tion and taking the limit of Ez <SC kT which eliminates 
the Boltzmann factors, the optical Bloch equations read 
using the basis states introduced in 



FIG. 12: (a) Example gatesweep. This plot has been obtained 
by increasing the gate voltage step by step and for each step 
taking a single-QD photoluminescence spectrum. The three 
strongest emission lines are identified as X", and ten- 

tatively X^"*" which result from s-shell electron-hole recom- 
bination from difi'erently charged excitonic complexes. Dis- 
crete and characteristic PL energies, different for every exci- 
ton complex, are induced by Coulomb interaction. Small con- 
tinuous PL energy shift is due to quantum confined stark shift, 
(b) Example differential transmission laser scan. When a laser 
is scanned over a QD transition light is resonantly Rayleigh- 
scattered. In DT experiment this Rayleigh-scattered light 
interferes with the laser background resulting in a dip in the 
transmitted light. The FWHM of the Lorentzian fit indicated 
by the solid red line is 460MHz. 

All experiments described here are carried out with a 
confocal microscopy setup immersed in a liquid helium 
bath cryostat at a temperature of 4.2 Kelvin. The nu- 
merical aperture of the microscope is 0.68 resulting in a 
diffraction limited spot size of ~ l^m. Area density of 
QDs in our sample is low enough to have < 3 dots in 
the focal spot simultaneously. Different QDs can then 
easily be spectrally separated by their inhomogeneous 
broadening. A magnetic field of up to 10 Tesla along 
the z-axis can be applied by a superconducting magnet. 
Piezo-electric nanopositioners allow us to move the sam- 
ple in XYZ-space. Transmitted light is collected and sent 
to a circular polarization analyzer which distributes the 
light to two photodetectors, similar to that of Ref.jl]. 

The initial step of our experiment is a gatesweep i.e. 
a PL measurement as a function of gate voltage. For 
this we send in a laser exciting electrons and holes in the 
bulk GaAs at an energy of «1.6 eV. After filtering out 
the pump light the resulting QD luminescence is sent to 
a grating spectrometer with a resolution ^ 20/xeV and 
detected by a hquid-nitrogen cooled CCD. 

Hereafter the differently charged excitonic complexes 
can be identified by their characteristic emission energy 
and voltage dependence profile [joj. From then on we 
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We have p'^^ = p[%, p'^^ = p'*^, p'^^ = P23' and pn + P22 
P33 = 1- 



APPENDIX C: DRESSED-STATE 
TRANSFORMATION 

The transformation used to diagonalize the coupling to 
the quasi-static nuclear (Overhauser) field can be written 
as 



H = SHS'' 
p = SpS'' 



(CI) 
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with (j) = ^h/^z and S"^ S = I. We assume (p 1. 
The new basis p3p determines the transformation ma- 
trix. When taking only first order terms 



S ■ 



The spontaneous emission terms then yield 




(C2) 



^ -{2Sa.2^S^ p Sa^^S^ (C3) 
For the new projection operator Sa23S^ we obtain 



(C4) 



and the conjugate relation. Here, tr^j = Using this 

with the previous relation, we obtain 



'S'irolaxation.rS'^ = - (2o-j3 p (Jg J — (T33P — pCTgg XC5) 

f . ^ ^ ^ , 

r _ 

-2^2 (^^23^^31 + <^i3P^32) 

where 7 = (jpT and F = F. aij = \i){j\ is the projection 
operator acting on p. The first, 7 term corresponds to 
relaxation via a weak optical transition induced by the 
hyperfine field, allowing for spin-flip Raman events and 
the second, F term, describes relaxation via the strong 
optical trion transition. 

The last term describes coherence induced by the spon- 
taneous relaxation into a superposition of dressed-basis 
ground states at a rate proportional to the occupation of 
the excited state ^33. When multiplying with (2| from 
the left and |1) from the right we obtain 



The transformed k terms keep the Lindblad form and 
we obtain for Ez e ^ kT 



relaxation, 



i2N2ipNi2 - Miip- pMii) (C7) 
+ {2Ni2pN2i - M22P - PM22) 



with N21 = Sa2iS\Ni2 = Nl^,Mu = Sa^S^ and k 
APPENDIX D: NUMERICAL STUDIES 



The derived formalism considers a static randomly ori- 
ented nuclear field. Within a measurement time -Bat 
changes ~100 times. In order to calculate measurable 
quantities such as linewidths and peak heights as func- 
tions of electric and magnetic field we have thus per- 
formed numerical simulations: For a given set of pa- 
rameters the steady-state solutions of the optical Bloch 
equations as given in Appendix |B] are numerically evalu- 
ated, in particular Im(p23(oo)) is then linked to the ab- 
sorption (details in section HV)) . A fluctuating hyperfine 
field is implemented by pulling three random numbers 
BN,i following (O. From Bjy^xy, using pT|) . state-mixing 
strength fin (|12p and pure Zeeman splitting Uz are calcu- 
lated before evaluating the density matrix steady-state. 
This procedure is repeated in order to average over '--^lOO 
random settings of the hyperfine field. In the cases the 
simulation could not be performed throughout the whole 
parameter space, we confirmed in key regimes that re- 
sults agree well with that of a static Overhauser field 
with equal magnitude in x,y,z: From ^ and l|12p we 
obtain for the rms- value of flnit) 



CjePB \ ^ 



3^ 
nuc 



(Dl) 



Here, the assumption for the observed absorption Q to 
be made is 



e B 



N,xy 
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( {B%,.yit)) 
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dt 



P21 



-2cj>rp: 



'33 



(C6) 



The same relation is obtained for p2i when multiplying 
with (1| and |2) respectively. 



i.e., the averaging over the absorption strength for differ- 
ent settings of the hyperfine field approximately equals 
the strength of absorption for the average field magni- 
tude, equal in x, y, z. 
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[50]feg/A^ ^ li^fj:^{A,f\^{R,)\^lf) 

jj^{A^ J(/ + l)) /{gefJ-s)'^', this result is obtained 
by assuming a box-like confinement potential, fur- 
ther transforming the sum into an integral and using 
the localization volume of the QD electron given by 
V-'=Jdr\i,{R,)f = 8/iNuo). 

[51] In this context we note that other experiments suggest 
that the dynamics of the nuclear spin ensemble could 
be altered by the measurement itself and the hyperfine 
field stays locked for a time on the order of seconds 0] 
Fig. 4(b) . In our measurements there is some evidence for 
alteration of dynamics at large magnetic fields and large 
gate voltage detunings (not shown in this work) ; however 
at low magnetic fields these locking effects do not seem 
to play a dominant role. 

[52] The type of co-tunneling relevant for us is inelastic co- 
tunneling as it leads to a transition between the spin 
ground states. 

[53] With Zeeman splitting 7^ ^^i^T i-^- the rate flipping 

the spin down is different from the rate flipping it up. The 
term /(e)[l — /(e)] has to be replaced by f{e^huJz/2)[l — 
fie ± hujz/2)\ for Kj,->t and K|->i respectively. 

[54] From the measured light polarization in our experiment 
we obtain, that the unwanted (T~-component is sup- 
pressed by a 25 X factor. An additional factor comes from 
the optical detuning with respect to the diagonal transi- 
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tion in magnetic field, with the laser being on the strong 
transition. At B =60mTesla, the detuning leads to a fur- 
ther suppression ~ 1/50, altogether suppressing the un- 
wanted excitation of the weak transition by more than 
1200 X, increasing with magnetic field. At fields less or 
on the order of the hyperfine field, spin ground states are 
strongly mixed leading to all possible (vertical as well as 
diagonal) couplings in the four-level scheme, thus creat- 
ing two differently polarized A-systems with equal decay 
rates to both ground states. The presence of the (un- 
wanted and mainly suppressed) a~ polarized light would 
just increase the bidirectional OSP by a small amount. 
Including Zeeman splitting, relation (I17|l can be writ- 
ten in the form = exp{ — fkJz/{2kT)) J^deA{e)f{e) 
and similarly — exp(fi,tJz/(2feT)) deA(e)/(e) 
with /(e) = [(1 + exp(e + hu^/ {2kT))){l + exp(e - 
hujz/{2kT)))]^^ . In (|32p . the term has to be replaced 
by whereas the K{n + 1) is replaced by 
As a consequence, the ratio of spin-up versus spin-down 
state occupation is governed, as one would expect, by the 
Boltzmann factor = exp{ — hu!z/{kT) 
This can be seen from (til^hmix |-ffint,rad| i) / as op- 
posed to iUit |H'i„t,rad| i) = 

Due to our square-wave modulation scheme inducing 
a non-exponential spin decay the conversion from the 



known modulation frequency to Kcxp rate does not need 
to include a factor of 2n. This has been confirmed by com- 
paring the result of a numerical simulation using square 
wave-shaped spin relaxation to the outcome of a simple 
rate equation model including exponential decay of spin. 
[58] As a Gedanken experiment we assume transfer of po- 
larization from light to nuclei via the electron spin un- 
der a static positive magnetic field. Then, due to nega- 
tive electron g-factor the spin-down state has higher en- 
ergy than spin-up state. The laser is resonant with the 
red (first case) or the blue (second case) Zeeman tran- 
sition. In the first case electron spin is pumped to the 
spin-up state which leads to nuclear spin-up polarization 
after an electron-nuclei fiip-fiop interaction. As the hy- 
perfine constant A is positive the energy of the electron 
spin-down state is consequently lowered while the trion 
states remain untouched, leading to a blue-shift of the 
red Zeeman optical transition which is being observed. 
The same happens in the second case where spin-down 
nuclear polarization is created and the energy of the spin- 
up state is lowered, again leading to a blue-shift of the 
observed transition. This scenario is clearly opposite of 
what we observe, thus ruling out the presence of an effi- 
cient DNSP. 



